Frustrated Ising model on the Cairo pentagonal lattice 
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Through the direct decoration transformation approach, we obtain a general solution for the pen- 
tagonal Ising model, showing its equivalence to the isotropic free-fermion eight-vertex model. We 
study the ground-state phase diagram, in which one ferromagnetic (FM) state, one ferrimagnetic 
(FIM) state, and one frustrated state are found. Using the exact solution of the pentagonal Ising 
model, we discuss the finite-temperature phase diagrams and find a phase transition between the 
FIM state and the disordered state as well as a phase transition between the disordered state and 
the FM state. We also discuss some additional remarkable properties of the model, such as the 
magnetization, entropy, and specific heat, at finite temperature and at its low-temperature asymp- 
totic limit. Because of the influence of the second-order phase transition between the frustrated and 
ferromagnetic phases, we obtain surprisingly low values of the entropy and the specific heat until 
the critical temperature is reached. 



I. INTRODUCTION 

Over the past six decades, much effort has been de- 
voted to determining the critical behavior of statistical 
properties of lattice models, which would allow a deeper 
understanding of order-disorder phenomena in magnetic 
solids. Following Onsager's pioneering exact solution for 
the square lattice Ising model exact solutions were 
also obtained for other regular two-dimensional lattice 
structures Q. In particular, exact results have been 
attained for the triangular, honeycomb, kagome, and 
bathroom-tile lattices [j], as well as for two- 
dimensional models, such as the Union Jack (centered 
square) Q and the square kagome 0] lattices. 

Geometrical frustration is mainly based on the trian- 
gle and tetrahedron structures, but it was also found in 
the Ising model on a pentagonal Penrose lattice proposed 
by Waldor et al.^ and solved exactly using the transfer 
matrix approach. 

More recently Urumov Q considered the Ising model 
on the Cairo pentagonal lattice using the decoration 
transformation [To| . This model has been mapped onto 
the Union Jack lattice [|| and its critical temperature 
and spontaneous magnetization properties have been dis- 
cussed. This model is interesting from a mathematical 
point of view. A few years ago, real materials with a 
Cairo pentagonal lattice structure were found; for exam- 
ple, the Fe 3+ lattice in Bi2Fe40g (described as a pen- 
tagonal Heisenberg model) was discussed by Ressourche 
et al. This material shows magnetic frustration. 

Also, theoretical calculations of the phonon structure 
of antiferromagnetic Bi2Fe4C>9 (space group Pbnm No. 
55, T « 240K) were studied using lattice dynamics and 
these results were confirmed experimentally by polarized 
Raman spectroscopy from 10 to 300 K [12|. More re- 
cently, some additional experimental studies were per- 
formed [l3|, Q3] ■ Ralko [lj| also discussed the hard-core 
extended boson Hubbard model on the Cairo pentagonal 
lattice, using the numerical quantum Monte Carlo study 



of stochastic series expansion and cluster mean-field the- 
ory. 

The purpose of this paper is to present a general exact 
solution of the pentagonal Ising model and as a special 
case we obtain the Urumov solution using a standard 
decoration transformation approach [10]. Furthermore, 
we present a more simplified solution through the direct 
decoration transformation [l6| instead of the standard 
one [l(|- The generalized version of the latter (l7l . [iH | 
is widely used to solve some two-dimensional decorated 
Ising [l9j and Ising-Heisenberg (20l - [22| models. 

This paper is organized as follows. In Sec. II we 
consider the detailed description of the Ising model on 
the Cairo pentagonal lattice. In Sec. Ill we discuss 
its phase diagram at zero temperature. Section IV is 
devoted to the pentagonal Ising model mapping, using 
the direct decoration transformation flH | for the isotropic 
free-fermion vertex model [23j . presenting the most rel- 
evant results and discussion. In Sec. V we obtain the 
finite-temperature phase diagrams, critical temperature, 
magnetization, entropy, and specific heat. Section VI 
summarizes our discussion. 



II. ISING MODEL ON THE CAIRO 
PENTAGONAL LATTICE 



The highly anisotropic Heisenberg model on the Cairo 
pentagonal lattice considered by Ressourche et al. (TlT | 
could be reduced to the Ising model on a Cairo pentag- 
onal lattice. Therefore, let us consider the Ising model 
on a planar lattice where the tiling is achieved with non 
regular pentagons; the lattice may be viewed as an as- 
sembly of checkerboard ordering with the elementary cell 
(see Fig. [5]) rotated by tt/2 in the neighboring square pla- 
quettes, as shown in Fig. Q](for more details see Ref. |0|). 

The Hamiltonian of the Cairo pentagonal Ising model, 
(represented schematically in Fig. [T]), discussed previ- 
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Figure 1: 
lattice. 



Schematic representation of the Cairo pentagonal 



ously by Urumov Q , is expressed by 



H 



(i,3) 



J £ S k Tl, 
(k,l) 



(1) 



where the first summation is the contribution of the inter- 
action between the nearest neighbor with spin Sj (sj inter- 
acting with coordination number 3) and J\ corresponds 
to the interaction between Si and Sj. While the second 
summation is the contribution of the nearest-neighbor 
interaction J between spin s k and spin t; (r/'s interact- 
ing with coordination number 4) , conveniently we assume 
,s ; = ±1 and t i = ±1. 



III. PHASE DIAGRAM AT ZERO 
TEMPERATURE 

In this section we discuss the phase diagram at zero 
temperature of the Hamiltonian given in Eq. ([1} . 

In order to discuss the phase diagram at zero temper- 
ature, we define the magnetization M for the pentagonal 
lattice that is used throughout the paper, given by 



M = 



M + 2A/i 



(2) 



with Mo = (ti) and Mi = (si). 

The energy per plaquette of three ground states that 
appear for the pentagonal lattice Ising model is expressed 
in terms of an elementary cell (see Fig. [5]). It is worth 
highlighting that the elementary cell should not be con- 
fused with the unit cell of the pentagonal lattice. 

(i) The ferromagnetic (FM) state or saturated state has 
a total magnetization M = 1 and ground-state energy 
per plaquette E = — J\ — 4 J. Thus the FM state can be 
represented as 
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Figure 2: Unit cells Ai (i — 1,2,3) with magnetization M 
1/6 and unit cell Si with M = -1/6. 



This state is limited by J > for J\ > and J\ > — J 
for J i < 0, as displayed in Fig. 0J 

(ii) The ferrimagnetic state (FIM) has a total magne- 
tization M = 1/3 and ground-state energy per plaquette 
E = —J\+ 4 J, which corresponds to the configuration 
displayed in Fig. [5] Analogous to the previous case, we 
describe the state by 



| FIM) 



+_ + 
+~ + 



(4) 



This state is limited by J < for J\ > and J\> J for 
Ji < 0, as illustrated in Fig. 0J 

(iii) The frustrated state (FRU) is given as a combi- 
nation of states | ) with its rotated elementary cell 
and spin inversion on the elementary cell with ground- 
state energy per plaquette E = J% — 2| J|, which can be 
expressed by 



| FRU) = combinations of { 



)}• (5) 



states 



This state is limited by J\ ^ — | J| (see Fig. g]). 

We define m as the magnetization for each frustrated- 
state configuration in the range between m = — 1/6 and 
1/6 ( we denote by M the average of total magnetiza- 
tion). Combining the state of the elementary cell dis- 
played in Eq. ([5]) and its rotation in 7r/2 of the ele- 
mentary cell, it is possible to generate the geometrically 
frustrated state. In particular, when we combine half 

^ and the remaining states with _+ + \, we 

obtain an antiferromagnetic state, with null total mag- 
netization. Other intermediate states with magnetiza- 
tion ^ m 1/6 also could be obtained by combining 
the elementary cell state with different relative amounts 
of the state given by Eq. ([5])- More specifically, the 
unit cell magnetization could be classified as displayed in 
Fig. [21 i.e, Aj (i = 1,2,3) and B\ 1 with magnetization 
m = 1/6 and —1/6, respectively. Certainly, this is not 
the only way to classify the unit cell by its magnetiza- 
tion; any other classifications of the unit cell lead us to 
the same kind of lattice configuration. In Fig. |3]we show 
two particular situations of such a lattice configuration, 
formed by the unit cells of type Ai and B. In the bulk 
limit we have a lattice with total magnetization 1/6. In 
Fig. (2Kb) the lattice is composed of unit cells Ai and B\ 
with a different concentration, generating a different to- 
tal magnetization of the lattice. Thus we obtain the total 
magnetization of each particular configuration. However, 
the average magnetization of all these configurations with 
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Figure 3: (a) In the top (left) panel, we illustrate schemati- 
cally an arbitrary lattice configuration, which is included in 
a frustrated state; the dashed square is shown explicitly in 
the panel below. ( b) In the top (right) panel, we assem- 
ble a lattice with two different unit cells A(m = 1/6) and 
B(m — — 1/6). In the panel below, the sublattice (the dashed 
square) is explicitly shown; the contributions of unit cells A 
and B are in different concentrations. 



equal energy will be null M = of the frustrated state 
in the interval — 1/6 < m < 1/6. 




Figure 4: Phase diagram at zero temperature of the pentag- 
onal Ising model. 



IV. THE PENTAGONAL ISING MODEL 
MAPPING 

Through the direct decoration transformation [l6| . it 
is shown that the pentagonal Ising model is equivalent to 



a Union Jack lattice, which is in turn mapped onto the 
isotropic free-fermion eight-vertex model [2J, [25| . From 
this equivalence we can obtain the thermodynamics of 
the pentagonal Ising model. In order to study the spin- 
1/2 Ising model on a pentagonal lattice, we introduce the 
notation for each elementary square plaquette depicted 
in Fig. [5] (this should not be confused with the unit cell). 
In this case we assume {ti,T2,T3,T4}=±1, {si, s 2 }=±1 
and a = ±1. 
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Figure 5: Schematic representation of the elementary plaque- 
tte and its mapping to the Union Jack lattice . 

The Hamiltonian for an elementary plaquette as dis- 
played in Fig. [5] is given by 

- ftH e = J[sis 2 + J'is^n + t 2 ) + s 2 (r 3 + n)}, (6) 

where we are assuming J 1 = /3Ji and J = ft J and ft is 
defined by ft — l/fc^T, with ks the Boltzmann constant 
and T the absolute temperature. 

Instead of applying the standard decoration transfor- 
mation [n| [T3j as applied by Urumov [j| , we apply the 
direct decoration transformation [1611 to transform the 
plaquette into the Union Jack lattice (6j. In this case the 
Hamiltonian associated with Union Jack lattice can be 
expressed by 

-ftH e =K' +K[a(T 1 +T2+T S +T A )+K' 2 (riT 2 +T 3 Ti). (7) 

The Boltzmann weights of the Hamiltonian ^ may be 
written 



>({t})= E eJ(Sl 



S2+J'[si(t 1 +T2) + S 2 (t3+T4)] 



(8) 



si ,s 2 =±l 



whereas the Boltzmann weights for the transformed pla- 
quette is given by the relation 

^({r}) _ c K;- l +K' 1 a{T 1 +T 2 +T3+T 4 )+K^{r 1 T2+T 3 r i ) _ ^ 

<7=±1 



Similar to the previous notation, here we consider K a = 
ftKo, K l = ftKi and K 2 = ftK 2 , where K Q is taken as 
a constant shift energy, K 1 is the interaction parameter 
between the internal spin a and each spin r, and K 2 is 
the coupling term between spins r. 

Using the direct decoration transformation proposed 
in Ref. [lf|, we need to impose the condition, w({t}) = 
w({t}) between Eqs. ((HJ and © for arbitrary r. There- 
fore, we obtain only four non equivalent configurations 
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{ti,t 2 , t 3 , r 4 } ={+, +, +, +}, { + , +, +, -}, { + , +, -, -}, 
and {+,—,+,—}; any other permutation or spin inver- 
sion falls onto one of these configurations. Thus the 
Boltzmann weight for each configuration is given by 

6 = w{+, +, +, +) =2c A 'o +2A ^ cosh(4/^), (10) 
t,2=w{+ ,-,+,-) =2e K '°- 2K '*, (11) 

(12) 



&=w(+,+,-,-)=2e K ° +3K ', 
6 = w(+, +, +, -) =2e K ° cosh(2i^), 



(13) 

where 6 = £4 and £5 = 6 = £7 = 6- 

The above equations satisfy the isotropic free-fermion 
condition [23] W1W2 + W3W4 = w^wq + w?w%, following 
the eight- vertex model with Boltzmann weights wi, ...,u>8 
(w should not be confused with w) displayed in Fig. [5] 
Hence the free-fermion condition may be rewritten in 
terms of £ as 



2£f = (£i + &)&. 



(14) 



Therefore, the Boltzmann factor of an effective Union 
Jack lattice can be expressed in terms of the pentagonal 
Ising model coupling parameters 



£l =ru 
6 =2 (r + r 



6 =2r + r- 1 u- A 



2r- 



6 



(u 2 + u- 2 ) 



(15) 
(16) 
(17) 
(18) 
J' 



For simplicity we used the notation r = e Jl and u = c" 

Due to the step by step decoration transformation per- 
formed by Urumov Q (see Fig. 2), J\ was restricted 
only to Ji > or r > 1 (ferromagnetic); otherwise, if 
we consider J± < 0, we get a imaginary parameter in the 
intermediate transformation proposed by Urumov. How- 
ever, using the direct decoration transformation, we do 
not have such a restriction (Fig. [5]); but only r > 0, which 
means that J\ exchange coupling could be ferromagnetic 
or antiferromagnetic. 

Hence the pentagonal Ising model is completely equiv- 
alent to the Ising model on the Union Jack lattice [f| 
with the isotropic nearest-neighbor interactions defined 
by K\ and non crossing diagonal interactions between 
the second nearest neighbor given by Ki- In contrast, 
the Union Jack lattice was mapped onto the isotropic 
free-fermion eight-vertex model [23] by Choy and Bax- 
ter [U. Therefore, we relate the Boltzmann factor given 
by Eqs. (jTCTj) - (| 1 3 p and the Boltzmann factor of the Union 
Jack lattice given by Eq. (4)ofRef. (24|. These relations 
are given by 



2£i 



0J 2 



26 



UJ 3 



2 and W5 = 



2£ 5 



(19) 

The schematic representation of the eight-vertex model 
is given in Fig. [6] 

These Boltzmann weights will be used in the following 
section to study the critical temperature and spontaneous 
magnetization of the Cairo pentagonal Ising model. 



Figure 6: (Color online) Eight-vertex model diagrams 



V. THERMODYNAMICS OF PENTAGONAL 
LATTICE 

In this section we discuss thermodynamical properties, 
such as the entropy, specific heat, and magnetization, as 
a function of temperature, as well as the critical temper- 
ature behavior. The thermodynamics of the pentagonal 
Ising model can be expressed following the results given 
by Fan and Wu [23] . The exact result for the free energy 
of the pentagonal Ising model is then given by 



1 

4tt 



2ji 



hi 



A{cf>) + y/QWj] <ty, (20) 



where 



A{4>) =\ {& + & + 266) + (6 - 6) v / 66 cos(0), 



Q(4>)= (6-6)v / 66cos(</») + i(6+6r 



66 466 -(6+6)' 



(21) 



(22) 



Once the free energy is known, we can obtain straightfor- 
wardly the critical temperature, magnetization, entropy, 
and specific heat. 



A. Critical temperature 

In order to study the spontaneous magnetization, fol- 
lowing the result obtained by Choy and Baxter [24| and 
using Eq. (|19p . the magnetization Mo = (tx) for spins 
with coordination number 4 is described by 



where 



M = v/l-P, 



26(6 + 6) 



6 



6 2 



266 



(23) 



(24) 



Equation (f2~4"|) is expressed in terms of the pentagonal 
Ising model Boltzmann factor. It is important to note 
that this relation is valid for arbitrary spins s± and s?, as 
shown in Fig. [5] Using the results obtained in Ref. [2J|, 
the critical point of the Union Jack lattice is obtained 
from the condition Wi — W2 = 2u>3 for w\ > W2 or W2 — 
w\ = 21U3 for W2 > Wi] in terms of k, this means that 
the critical points occur at k = 1. Equivalently, using 
the pentagonal Ising model Boltzmann factor £, we have 



(6 - 6) = 466 



(25) 
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This condition must satisfy the critical point. Rewriting 
Eq. (|25p in terms of r and u, the critical points must 
satisfy the relation 



r V2u, ^ + 2,g + (l-^)V2 

where r c and u c denote r and u evaluated at the critical 
temperature T c . The same expression could be obtained 
fromEq.(8) oftheRef. @; however, due to the standard 
decoration transformation [l(|[I3 used by UrumovQ for 
this model, we have to eliminate the intermediate param- 
eter Q. Once the intermediate parameter Q is eliminated 
from Eq. (8) of Ref. Q, it becomes identical to our 
Eq. (|2"rJ)) for the case of J\ > (ferromagnetic coupling), 
which was previously studied by Urumov 

In Eq. (J2H]) we provide a closed expression for the 
critical point of the pentagonal Ising model. The curves 
where the critical points occur are illustrated in Fig. [JJ 
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Figure 7: Critical points curve fo the pentagonal Ising model 
as a function of the parameters J/T c and Ji/T c . 
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Figure 8: Phase diagrams as a function of temperature T c and 
the parameter \J\ for two different values of Ji. In (b) the 
dashed lines show the limiting values of J — > and T c — > 0. 



The finite-temperature properties of the system arer 
investigated by considering the effect of parameters J and 
J\ on the critical behavior. In Fig. [JJwe display the crit- 
ical point regions in units of critical temperature, where 
the phase transitions of the FIM region to a disordered 
phase (DP) and a DP to a FM phases are illustrated. It 
is important to highlight that Fig. [JJ becomes similar to 
Fig. |1 when T c ->• 0. 

An alternative phase diagram is depicted in Fig. [SJ 
where the phase diagrams are illustrated in the (T c , \ J\) 
plane (since T c is invariant unde ther J — > — J exchange) 
for fixed parameter J\. In Fig. [8ja) the second-order 
phase transition line is shown in the (T c , | J|) plane when 
the parameter Ji is fixed at J\ = —1.0; in this case there 
is a DP. Concretely, for J < we show two regions: the 



FIM phase and the DP; for J > we have a DP and a 
FM phase. 

In Fig. [5Jb) we show the behavior of the critical 
temperature when J\ = 1.0, displayed BjS cl solid line. 
For low values of T « 0, the left-hand side of Eq. 
(I26p goes to infinity; this implies that the denomina- 
tor on the right-hand side must satisfy the condition. 
ul 2 — bu & c — 5u^ + 1 = 0. Thus we obtain the solution 
T c = ±2.2691J (dashed lines), which is valid within the 
limit J — > and T c — > 0. This result is in agreement with 
the phase diagram at zero temperature (see, for instance, 
Fig. |4|); more specifically, for J\ = 1, the phase transition 
occurs at J = 0. 

We now comment on the finite-temperature phase dia- 
grams displayed in Fig. O in which the critical tempera- 
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Figure 9: Phase diagram in the (T c /J,Ji/J) plane for the 
pentagonal Ising model. 



ture T c / J is shown as a function of the parameter Ji/J. 
Using the equation of the critical points (|26p . we obtain 
the plot illustrated in Fig. [9l From that we can analyze 
three limiting cases. 

(i) For j Owe obtain T c /J = 1.3084 (T c /J = 
— 1.3084) and the pentagonal lattice is reduced to 
a ferromagnetic (ferrimagnetic, with total magne- 
tization M = 1/3) decorated square lattice. 



(ii) For 4 



oo, from our calculation, for J > and 
Ji > 0, we find the solution T c /J = 2.2691; in 
this case the pentagonal lattice is reduced to the 
ferromagnetic square lattice. In the case J < and 
J i < , in the limit under consideration (4 — > oo), 
the pentagonal lattice falls into the bottom righ- 
hand corner, which is a DP state (see Fig. |5J). 

(iii) For 4 -> -oo we obtain T c /J = -2.2691 for J < 
and Ji > and the pentagonal lattice is reduced 
to a ferrimagnetic (with total magnetization M = 
1/3) square lattice. Meanwhile, for ^j- — > — oo for 
J > and Jj < the pentagonal lattice falls into a 
disordered state, which corresponds to the top left 
corner Fig. [9l 



When K-2 = 0, according to a pentagonal lattice map- 
ping onto an effective square lattice (see Fig. [5]), within 
the limit ^j- — > oo, the effective lattice is reduced to a 
ferromagnetic square lattice. In the limit ^j- — > — oo the 
pentagonal Ising model reduces to an antiferromagnetic 
square lattice. 



Figure 10: (Color online) Internal energy U as a function of J 
for a fixed value J\ — —1.0. The dashed line corresponds to 
the DP ground energy and the dash-dotted line corresponds 
to the FM ground energy. 



B. Internal Energy 

The internal energy of the pentagonal lattice Ising 
model defined by U = T 2 d ^J^ can be obtained straight- 
forwardly from Eq. (|20D . In Fig. [TU] we display the inter- 
nal energy as a function of coupling parameter J in the 
low-temperature limit in order to observe the low-lying 
energy contribution close to the critical temperature as- 
suming 3\ = —1.0. The dashed line corresponds to the 
DP internal energy at T = (the ground-state energy), 
given by U = — 1 — 2 J, while the dash-dotted line indi- 
cates the internal energy U = 1—4 J from which the phase 
transition at zero temperature occurs at J = 1. In the 
case T c — > and assuming J = 1 + 5 with S > 0, we have 
u c — > oo and then it is possible to write u c = r^ 1- * 5 . Fur- 
ther, by substituting into Eq. (p?S)) , after some algebraic 
manipulations we obtain r 2S sa where r c = 



was defined previously. Thus a second-order phase tran- 
sition occurs at T c w . It is worth highlighting 

ln(2V2) ° ° ° 

that the lowest critical temperature occurs at T c = 0. 
Therefore, the contribution of the low-lying energy is ab- 
sorbed by the second-order phase transition. Certainly 
there is no second-order phase transition for J < 1 and 
J i = —1. The solid blue (thick) line and the dark-blue 
(thin) curve represent the internal energy in the disor- 
dered phase, while by the red (thick) line and the orange 
(thin) line curve represent the ferromagnetic region for 
two critical temperatures T c = 0.192 and 0.385. 

In Fig. [TT] we plot the internal energy as a func- 
tion of temperature for several values of J around the 
second- order phase transition, assuming a fixed value for 
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Figure 11: (Color online) Internal energy U as a function of 
T and a fixed value Ji = —1.0. The dotted curve corresponds 
to U(T C ). 



Ji = —1.0. The black dotted line represents the inter- 
nal energy U (T c ) evaluated at the critical temperature T c 
given by Eq. (f5T|) at low temperature. Below the critical 
temperature, the internal energy is almost constant (fer- 
romagnetic phase), which means that it is mainly given 
by the zero-temperature ground-state energy. Although 
at a critical temperature there is a sudden change of cur- 
vature, this change becomes dramatic for lower critical 
temperature; for higher critical temperature this change 
of curvature becomes smoother. The lowest critical tem- 
perature occurs at T c = for J = 1; therefore, for lower 
values of the coupling parameter J there is no second- 
order phase transition. For a sufficiently high temper- 
ature the internal energy leads to an asymptotic limit, 
whereas the internal energy increases almost proportion- 
ally to the temperature. 



C. The spontaneous magnetization 

The total magnetization M is given by Eq. @ for the 
pentagonal Ising model, which will be discussed in order 
to show the spontaneous magnetization. It ought to be 
pointed out that the calculation of the magnetization of 
internal spin Sj, i = 1, 2, may be obtained following the 
results obtained by Choy and Baxter (24|, which could 
be expressed by the relation 



where after some algebraic manipulations the coefficients 
become 



ci 



(-2 



1 /Ei 

4 Ui 
1 

4 



— - 2— 

6 6 



(28) 
(29) 



Defining Si and £5 in analogy to the Boltzmann factors 
£, we have 



Si 

£5 



-2r 



-2r 



2ru 4 



2ru 



(30) 
(31) 



In order to obtain the three-spin correlation function 
(T1T2T3) we use a checkerboard Ising model equivalent to 
that used by Choy and Baxter (24|. Here we use the rela- 
tion obtained in Ref. (24|, but in our case we rewrite this 
relation in terms of the pentagonal Ising model Boltz- 
mann factor. Hence, using some algebraic manipulations 
we have 



(tit 2 t 3 ) = R(r,u)(n), 



where 
R(r, u) 



2£i 



fi-& 



£l±|3 I 1 _ 



(32) 



(33) 



Therefore, the magnetization M\ can be expressed as a 
function of magnetization Mq, which is given by 



1 



Mi = - 



Mr, 



Mi 



si 



Cl(T) + C 2 {TlT 2 T 3 ) 



(27) 



£1 £5 

(34) 

It is important to highlight that Eq. (|34|1 is expressed in 
terms of the original parameters of the pentagonal Ising 
model instead of parameters of the effective Hamiltonian 
such as those obtained by Urumov Using Eq. (p4| . 
we are easily able to manipulate the parameters of the 
pentagonal Ising model in order to discuss the sponta- 
neous magnetization. At the critical point we need to 

substitute the expression r = r c T , where r c is defined 
in Eq. ([26]). From Eqs. (J23J) and (JHIJ we can obtain a 
closed expression for the total magnetization of the Ising 
model on a pentagonal lattice, using the relation given 
byEq. ©. 

We now discuss the behavior of the total magnetiza- 
tion of the pentagonal Ising model as a function of the 
parameter J for the low-temperature limit. In Fig. I12l we 
plot the magnetization at low temperature as a function 
of J, where we display two types of plateaus for the FIM 
state and the FM state. This is in agreement with the 
phase diagram displayed in Fig. 2] and [71 whereas the in- 
termediate state corresponds to the FRU phase and the 
DP, respectively. 

Hence, in Fig. [T^a). for J x = -1.0, when T = there 
are three well defined regions: FM phase with M = 1. 
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Figure 12: Total magnetization of the pentagonal Ising model 
for three different values of T as a function of the parameter 
J and fixed parameter Ji. (a) Ji = —1.0 and (b) J\ = 1.0. 



the FIM phase with M = 1/3, and the intermediate FRU 
phase; also, by increasing the temperature (for example, 
from T = 0.1 to 1.0) the disordered phase increases (the 
| J | increases). Meanwhile, in Fig. [T27b). for J\ = 1.0 and 
at zero temperature, we have a direct phase transition 
between the FM state and the FIM state. However, for 
a nearly zero temperature T = 0.1, a small intermediate 
region arises that corresponds to the disordered phase 
region. For higher temperature such as T = 1.0, the DP 
region is even larger. 

Another way to analyze the total magnetization is by 
exploring the temperature dependence T/T c of the total 
magnetization. In Fig. [T3] we plot the magnetization as 
a function of temperature for two values of J/T c = ±1 
and J\ = —1.0, where the magnetization of the saturated 
FM state and the FIM state is illustrated. For the case 
of J/T c = 1 we have the FM region with M = 1 the 
total magnetization at zero temperature. The total mag- 




Figure 13: Temperature dependence of the total magnetiza- 
tion of the pentagonal Ising model for two different values of 
J/T c and a fixed value Ji = —1.0. 



netization vanishes at T = T c as temperature increases; 
therefore, a second-order phase transition occurs at the 
critical temperature. 

For J/T c = —1, the total magnetization corresponds 
to the FIM region, with M = 1/3 at zero temperature, 
and it vanishes at T = T c as displayed in Fig. [T3] Thus a 
second-order phase transition occurs again at the critical 
temperature, which is in agreement with the critical point 
curve displayed in Fig. [7] 



D. Entropy 



The entropy can be easily obtained as a negative tem- 
perature derivative from the free energy (|20| S = — ^ , 
while the specific heat can be written as a temperature 
derivative from the entropy C = T^. In what follows 
we consider only the case J > since for J < we have 
the same behavior. 

In Fig. [T3] we display the entropy as a function of 
temperature for several values of J with J\ = —1.0 
fixed. Figure [14] (a) shows the low-temperature behav- 
ior of entropy, where the residual entropy appears at 
S a = ln(2)/2 = 0.3465 for |J| < 1.0. This means that 
we are in a geometrically frustrated region, which is in 
agreement with the illustration of the phase diagram in 
Fig. [U The residual entropy is proportional to In (2); 
this number comes from the two configurations given in 
Eq.lp)]). For \J\ = 1.0 the residual entropy has a different 
nontrivial value as displayed in Fig. I14f b) . To obtain the 
entropy explicitly, we return to Eqs. ([2lTl and |22|) and 
set r> = — = e^ 1 l T ". Thus we obtain 



V2cos((f>) + 0{r 2 c ) 



9 



and 



-7 = 0.1 7=0.2 7=0.5 



rl 2 Q 



- + V2cos(0) + 2 cos 2 (0) + 0{rl). 



Finally, for T c -> (r c -> 0) and using Eq. flU), the 
residual entropy becomes 

1 r 27r 

S = — ln(| + |V2cos(0) + i| + v / 2cos(</)))d0 
47r Jo 

0.5732714757. (35) 

This is due to the degeneration of the phase boundary 
between the FRU and FM (or FIM) regions at T = 
(for detail see Fig. 3]). This result was derived in a way 
similar to that discussed by Wannier j2(| for the case 
of a two-dimensional triangular lattice. While in Fig. 
rUT c) there is no residual entropy for \J\ > 1, the stan- 
dard temperature dependence of entropy appears with 
a strong change of curvature located at critical points 
where a second-order phase transition occurs. 

It ought to be highlighted that in the low-temperature 
limit (below the critical temperature) the entropy (in the 
FM state) for J > 1 and with J\ = — 1 can be obtained 
from Eq. (|2"0")) . More explicitly, by fixing J = 1 + 5, where 
S > 0, it is possible to write u = — , with r = c~ and 

s = c^r. Substituting Eqs. and (J2S]) into Eq. (gQJ) , 
the integration of Eq. (|20| results in 



-3 - 26 
« 1-4J-2 
Finally, the entropy 5 



T(-21ns + 2rV) 

2(2| J|-l) 



^ can be written as 



S : 



2 A -A/T 

T 



(36) 



where A = 2(2| J| — 1) is the energy gap. In Figs. [T4T c) 
and [n| d) the low-temperature limit is well fitted by the 
above limiting expression. 

An additional plot of entropy S against J in the low- 
temperature limit is displayed in Fig. 1151 where the 
residual entropy is illustrated by the dashed black line 
at zero temperature. For J < 1 there is residual en- 
tropy S = ln(2)/2, while for J = 1 the residual entropy 
becomes S = 0.5732714757; for higher values of J > 1 
there is no residual entropy. Thereafter, we observed the 
entropy in the low-temperature limit, where we can show 
the effects of residual entropy. The low-lying energy con- 
tribution for the entropy between the DP and the FM 
phase is absorbed by the second-order phase transition 
as a consequence of the entropy falling dramatically to 
zero entropy for J > 1; for higher temperature the en- 
tropy change curvature becomes softer. 



E. Specific heat 

Finally, we conclude our analysis of thermodynamics 
by exploring the temperature dependence of the specific 
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Figure 14: (Color online) Entropy as a function of tempera- 
ture for Ji = —1.0 and (a) and (b) J ^ 1 and (c) and (d) 
J > 1 
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Figure 15: (Color online) Entropy S as a function of the pa- 
rameter J for a fixed low temperature and Ji = —1.0. 



heat. Some typical thermal variations of the specific heat 
of the pentagonal Ising model are plotted in Fig. [16] for 
several values of J and Jj = —1.0 fixed. In Figs. [ToTa) 
and (b) we present the temperature dependence of the 
specific heat in the DP states and show that there is no 
phase transition at finite temperature because we are ob- 
serving the frustrated region. In Figs. ITUf c) and (d) the 
specific heats are logarithmically divergent at the criti- 
cal temperature, which is associated with a continuous 
phase transition between the spontaneously ordered and 
disordered phases. Clearly, this means that we are facing 
a FM region, which can be verified in the phase diagram 
illustrated in Fig. 0J 

Plotted in Fig. ICTc) when J > 1 and J\ = — 1 is the 
specific heat versus temperature; a surprisingly almost 
null specific heat is displayed, until the critical temper- 
ature is achieved. When the absolute value of the ex- 
change interaction is only slightly above 1, the critical 
temperature (an order-disorder phase transition) can be 
obtained easily from Eq. (f2T))k this critical temperature 
is approximately given by the expression 



2(| J| - 1) 
ln(2\/2) 



(37) 



in the low-temperature limit. 

In Fig. [T7]we display the specific heat as a function of J 
in the low-temperature limit, where we show the specific 
heat behavior around the second-order phase transition. 
As discussed previously, the low-lying energy contribu- 
tion is absorbed by a second-order phase transition in 
the case of J > 1, while for J < 1 there is no second- 
order phase transition; then we can observe the low-lying 
energy contribution as a small anomalous broad peak. 




-J=0.9 7=0.95 7=0.97 7= 1.0 



1 

0.8 
0.6 



0.4- 



0.2 




0.5 



1 

T 



1.5 



-7=1.0- 



-7=1.02- 



-7= 1.1 ■ 



-7=1.2 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.! 
T 



-7= 1.35 7= 1.6 



d) 




I 


/ 

/ 

/ 



0.5 1 

T 



1.5 



Figure 16: (Color online) Specific heat as a function of tem- 
perature for Ji = —1.0 and (a) and (b) J ^ 1 and (c) and (d) 
J > 1. 
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Figure 17: (Color online) Specific heat as a function of J in 
the low-temperature limit and Ji = —1.0. 



Similarly, the low-temperature asymptotic limit for 
specific heat can be derived from Eq. (|20[) so that 



J = 1.02 J= 1.1 J= 1.2 




T 



Figure 18: (Color online) Low-temperature limit specific heat 
against temperature for the same values as in Fig. I16f cl. 
The solid line corresponds to the exact specific heat and the 
dashed line is the low-temperature limit of the specific heat. 



C « f^c- A/T (38) 

The energy gap is large enough even when J = 1.02 
because the order-disorder transition occurs for |J| w 
l + ^fln(2\/2). 

In Fig. [15] we show the magnification of Fig. ITBT c) 
in low-temperature limit, which is well fitted by Eq. (|38p . 
The solid line corresponds to the exact specific heat and 
the dashed line represents the low-temperature asymp- 
totic limit of the specific heat. For J = 1.02 the cor- 
responding low-temperature approximation is valid for 
T < 0.0385, while for J = 1.1 and J = 1.2 clearly the 
low-temperature curve accompanies quite well the exact 
solution. 



VI. CONCLUSIONS 

Using the direct decoration transformation (l6| . we 
have solved the pentagonal Ising model with a more gen- 
eral coupling parameter and compared it with Urumov's 
Q solutions. We have found a frustrated phase of the 
pentagonal Ising model. 

In addition, we have obtained a simplified solution for 
the free energy, as well as a closed expression for the 
critical temperature. Although this model has already 
been solved by Urumov through the standard decora- 
tion transformation [Tol \vi\ in the nonfrustrated region 
(Ji > 0), such a result contains unnecessary interme- 
diate parameters that can be avoided so that a closed 
expression similar to Eq. (|2"6")l can be obtained. We have 



studied the ground-state phase diagram, which exhibits 
a ferromagnetic state, a ferrimagnetic state, and a frus- 
trated state at J\ < —\J\. Following the exact solution 
for the pentagonal Ising model, we have discussed the 
finite-temperature phase diagram, as shown in Figs. [8] 
and[9j identifying the phase transition between the FIM 
state and the DP state and also between the DP state 
and the FM state. 

The analysis of the limits in Fig. |9] allows one to find 
three relevant phases. For J\ — > the pentagonal lat- 
tice reduces to a ferromagnetic (ferrimagnetic) decorated 
square lattice. For J\ — s- oo and J > the pentago- 
nal lattice reduces to the ferromagnetic square lattice. 
Finally, for Ji — > — oo and J < the pentagonal Ising 
model reduces to a ferrimagnetic square lattice. The to- 
tal magnetization as a function of the parameter J and 
for a fixed value of J\ for the ferromagnetic state [M = 1) 
and the ferrimagnetic state (M = 1/3) is shown in Fig. 

EH 

For a fixed value of J\ = — 1 there is a residual entropy 
S = 0.3465. For |J| < 1.0 and |J| = 1.0 a nontrivial 
residual entropy Sq = 0.5732 is found, as shown in Fig. 
IT4l Because of the frustrated state, the entropy below 
the critical temperature shows a strong change of cur- 
vature for J > 1. The specific heat capacity was also 
investigated at fixed J\ = — 1 and |J| < 1 [see Figs. 
[TTjfa) andfTBTb)]. For J\ = —1 and | J\ > 1 we have un- 
usual behavior due to frustration of the entropy and the 
heat capacity at temperatures below the critical value, as 
shown in Figs. ITBT c) andfTHTd). 
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